A pseudo-potential analog for zero-range photoassociation and Feshbach resonance 
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A zero-range approach to atom-molecule coupling is developed in analogy to the Fermi-Huang 
pseudo-potential treatment of atom-atom interactions. It is shown by explicit comparison to an 
exactly-solvable finite-range model that replacing the molecular bound-state wavefunction with a 
regularized delta-function can reproduce the exact scattering amplitude in the long- wavelength limit. 
Using this approach we find an analytical solution to the two-channel Feshbach resonance problem 
for two atoms in a spherical harmonic trap. 
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Coupling between atoms and molecules in quantum- 
degenerate gases is an ever-present aspect of ultracold 
atomic physics. Feshbach resonances (FR) Q] are now 
routinely used for control over atomic interactions Q 
and the formation of molecular Bose-Einstein conden- 
sates (3- ED- 0^ E3l Laser- induced photoassociation (PA) is 
also widely employed 0, , having the advantage of con- 
trol over the coupling strength j9j. While a zero-range 
approach to atom-atom collisions has long been a cor- 
nerstone of BEC theory, an analog to the Fermi-Huang 
pseudo-potential approach 10] has yet to be formulated 
to treat multichannel free-bound coupling in ultra-cold 
atomic gases. 

In the long- wavelength limit, the energy-dependence 
of the scattering phase-shift for atomic collisions takes 
a universal form, with all information about the details 
of the interaction potential contained in a single param- 
eter, the scattering length. As a result, the full inter- 
action potential can be replaced by a regularized delta- 
function pseudo-potential, which yields the correct scat- 
tering amplitude up to a third-order correction in the 
ratio of the effective range to the incident wavelength. In 
this Letter we formulate an analogous approach to atom- 
molecule coupling, by replacing the bound-state wave- 
function with the zero-range object which correctly re- 
produces the long-wavelength scattering amplitude. The 
resulting model contains no divergences and does not re- 
quire a momentum cut-off. It is likely that this model will 
play an important role in understanding the role played 
by atom-atom correlations in FR and PA physics, partic- 
ularly in the strong-coupling regime, where such effects 
play a dominant role. 

The first zero-range model for BEC atom-molecule 
coupling, proposed by Heinzen and coworkers, replaced 
the bound-state wavefunction with a delta-function . 
This approach was shown by Holland and coworkers to 
contain a UV divergence when pair correlations were 
taken into account |l2j, thus limiting its applicability. 
Holland and coworkers demonstrated that this diver- 
gence could be removed via a momentum cut-off and 
re-normalized detuning. As we will see, this approach 
fails in the presence of a background scattering length. 

We begin our analysis by considering a pair of atoms 
described by a relative wavefunction <pj(r,t), where j = 



1, 2 corresponds to an internal spin state. The eigenstates 
of this system obey the Schrodinger equation, 



Ect> j (r) = - — V 2 <j> j (r) 



J2V jk (r,t)4, k (r), (1) 



where E is the energy eigenvalue, fi is the reduced 
mass and Vjfc(r) is the inter-atom potential. For our 
model system we assume that the first channel sees a 
flat potential, Vu(r) = 0. The second channel sees 
a spherical- well potential of depth Vo and radius w, 
V22(r) = Uo — VqU(w — r), where Uo is the continuum 
threshold energy and U(x) is the unit-step function. In 
the absence of coupling terms, i.e. for Vvz(r) = 0, the 
spectrum of the second channel consists of a continuum 
of states above the threshold energy, Uq, and a discrete 
set of bound states with energies between Uq and Uq — Vq. 
The bound-states are all of the form 



Mr) = 



: r > w 



u sm(kbW) r 



(2) 



where ab and kb satisfy the equations |p \k 2 + 1/a 2 ] = 
Vo and cot(fcbiu) = — l/(fe(,a(,), an d A/j, is determined by 
normalization. The bound state energies are Eb = Uo — 

We proceed by first expanding the second channel 
wavefunction, 2 (r), onto its bare eigenstates under the 
simplifying assumptions that only a single bound state 
is near-resonantly coupled to the first channel so that all 
other states may be neglected. We assume the interac- 
tion potential has the form V\2(r, t) ~ L-H e -*w*. Taking 

E = h 2 k 2 /(2ii) then leads to an eigenvalue problem for 
a continuum coupled to a single bound state, 



1 



[k 2 +V 2 ] <Mr) = GM?)c 



(3) 



~ [k 2 - 2A] c = G* j d 3 rr b (r)Mr), (4) 
where c is the probability amplitude for the atom pair 



to be in the bound state, and A 



-(U - to) 



l 



is the detuning away from the atom-molecule resonance 



2 



at k = 0. The coupling constant G will depend on the 
details of the atom-molecule coupling scheme. 

Our goal is now to solve this eigenvalue problem, under 
the boundary conditions 



lim 4>\{y) = 



p — ikr Akr 
~ +/ — 

r r 



limr0i(r) = 0, 

J >0 



(5) 
(6) 



in order to determine the scattering amplitude / = /(fc). 
The solution can be obtained via the ansatz 



(r) 







+ f e - 

sin(fcr) 



2Gaj 
l + (a b fc) 2 
2G 



+ k2 ^ K 2 Cipb(r) :r<w 



(7) 

This ansatz explicitly satisfies ©, as well as the bound- 
ary conditions (151011 . Equation (0J, together with the 
continuity equations <fii(w + ) — (fti(w~) and V0i(u> + ) = 
V4>i(w~) can then be used to determine the three un- 
knowns /, c, and [3. These equations are linear in the 
three unknowns, and can be thus solved in a straightfor- 
ward manner. 

The long- wavelength limit requires that 1/fc be large 
compared to the size of the bound-state. As the size 
of the bound-state is w + af,, this is equivalent to the 
limits kw <C 1 and kcib <C 1. For our model potential 
the condition Kb > l/w is always satisfied, so that k/Kb 
is a small parameter as well. Expanding the scattering 
amplitude f{k) in terms of these small parameters then 
yields 



f(k) 



25 - ik^£ 



k 2 



25 + ik 



0\ 



(8) 



where e € {kw,kab,k/Kb}, and we have introduced the 
light-shifted detuning 



A - 8ir\G\ 2 N?e- 2w/a 
1 



1 



K 2 a 2 



al(a b 



e 2w/a 



■ w) 



1 + 



K 2 b a 2 



(9) 



and the effective coupling constant 



X = 4irGNbe- w/a a b (a b + w) 



Kla 2 



(10) 



The important point here is that all of the details of 
the potential can be absorbed into effective detuning and 
coupling constants. 

We now consider a zero-range model in which the 
bound-state wavefunction V'f)( r ) i n © is replaced by a 
regularized delta-function, GV>b(r) — * x^{ r )~§p T - l n a d- 
dition, the detuning A is replaced by the light-shifted 
detuning S and the coupling constant G is replaced by 



the effective coupling constant x- The Schrodinger equa- 
tion for this model is given by 



i[fc 2 +V 2 ]^(r) = X S 3 (r)c 
\ [k 2 ~ 26] c 



(11) 



d 



X* I d 6 r6\v)— r^i(r). (12) 



This problem can be solved by making use of the ansazt 
x (r) = £-^1 + fi^- and the identity V 2 ± = -47n5 3 (r). 
The scattering amplitude is readily found to be 

k 2 -25-ik^- 



f(k) 



k 2 -25 + ik&£ 



(13) 



which agrees with the result (|SJ up to a correction 
of third-order in the small parameters kw, kab, and 
k/Kb- Thus the zero-range model <|11I12[1 will reproduce 
correctly the long-wavelength atom-molecule quantum- 
dynamics of our model potential. 

Second quantization of this model yields the Hamilto- 
nian 



n 



4m 



<Tr 



ft(r)V 2 i;{r) + -#t(r)V 2 tf(r) 



V2m 
+H.c. 



[ d 3 Rd 3 r&(R)5 3 (r)^riJ>CR+-U(-R- 
J or 2 



"-) 
2' 

(14) 



where ip(j) is the annihilation operator for an atom of 
mass m = 2/i, and ^(r) is the annihilation operator for 
a molecule of mass 2m. The system of equations <|1U 
1121 can be derived from this Hamiltonian via the 2-atom 
quantum state 



|*> = —= / d 3 Rd 3 r$(R)(t)(r)^(R+ -)^(R- 
v 2 J 2 



d 3 m(R)¥(R)\Q), 



5 J\0) 
(15) 



where $(R) is an arbitrary center-of-mass wavefunction 
and |0) is the vacuum state. The Hamiltonian 1|14|1 should 
form the basis of any field-theoretical description of zero- 
range atom-molecule coupling. 

As an example, we now solve the problem of two 
bosonic atoms in a spherical harmonic oscillator (with 
frequency LOtrap) with both s-wave collisions and cou- 
pling to a bound state in a second channel. With 
E = TiLj trav (y n + 3/2), 5 — > %ujt rap 5, and using harmonic 
oscillator units, the time-independent Schrodinger equa- 
tion can be written as 

\v n + ±V 2 - \r 2 + |] 0„(r) = 2n (f ) 5 3 (r)-§- r r 0„(r) 

+^ 3 /W(r)c„ (16) 
K -5}c n = n 3 / 4 n J d 3 r 5 3 (r)JLr 0„(r), (17) 

where n is an integer label for each quantum level (the 
lowest energy level corresponding to n = 0), a is the 



3 



background scattering length, A is the harmonic oscillator 



length of the trap, and SI = A 2 7T 
eigenfunctions are found to be ^ 



3 / 4 X- The normalized 



<t>n(r) 



27r3/^K,f) < 



si 2 ^ v[-if] [^{-^)-i>{- 



2 



)] 



r[-^ti] 



(18) 



-1/2 



where = 1 - 2aT[-f]/r[ 



2 



(19) 
U(a,b,z) is 



the confluent hypergeometric function and i}>(z) is the 
polygamma function [14j. The eigenvalues {^ n } are de- 
termined by the characteristic equation 



r[-i^]^ n a)' 



(20) 



where there is an apparently non-trivial relation |c„| 2 = 
du n /dS. It is straightforward to show that the spectrum 
of eigenvalues will agree exactly with those of a single- 
channel system with the energy-dependent effective scat- 
tering length 



eff{v) = a + o 



2{ v -sy 



(21) 



which is the familiar Feshbach Resonance result. The 
only difference between the true atom-molecule eigen- 
states and the equivalent single-channel states with scat- 
tering length a e ff, is the presence of the bare- molecule 
population, |c„| 2 . From a series expansion of (|18|) the 
l/r part of (j) n (r) is found to be - 27r i/ 4 "(i/„,") 0nly 
for a — is this term independent of v n , so that it can 
be removed via a renormalized detuning |12j . 

On resonance we have v n = 5 and |a e //| — > oo. A 
careful analysis shows that this requires ^„ = 2n — 1 
and c„ 7^ 0. Thus the eigenvalues are driven to odd- 
integer or 'fermionized' values, for which the regular part 
of 0n( r ) vanishes at r = 0. Inserting this result into Eq. 
(|19fl gives an analytic expression for the on-resonance 
molecular fraction, 



1 



i + a„n 2 '■ 



(22) 



where a n 



T 7r/2. For the low lying levels we have 



(2n-i)!! 

aa = 7r/2, ol\ = 7r and ot2 — 47r/3. 

The energy-dependence in the effective scattering 
length is critical to understanding the cross-over between 
the weak-coupling and strong coupling regimes. The re- 
quirement for a significant deviation from the bare-trap 
spectrum is o e ///A ~ 1. Obtaining this condition via 
Feshbach resonance requires S = v n ± y^fi 2 /2. If this 
width is smaller than the level spacing, only a single level 
can be near-resonant for a given detuning. In this weak- 
coupling regime, il 2 <C 1, the spectrum consists of a series 




FIG. 1: The eigenvalue spectrum as a function of the detun- 
ing for the case a = .3A and SI = .2, illustrating a sequence of 
avoided crossings in the weak-coupling regime. The dashed 
lines correspond to the uncoupled eigenvalues 
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FIG. 2: The eigenvalue spectrum as a function of the detun- 
ing for the case a — and f2 = 10, illustrating the 'fermion- 
ization' of the low lying levels in strong-coupling regime (in 
the vicinity of S = 0). The dashed lines correspond to the 
odd-integer values v n =2n — 1. 



of avoided crossings between the bare molecular level and 
the uncoupled eigenstates of the 'open' channel. At each 
avoided crossing there will be strong mixing between a 
single trap level and the molecular state. Sweeping the 
detuning can select which trap level is resonantly cou- 
pled to the molecular state. This is illustrated in Figure 
n where we have plotted the eigenvalue spectrum as a 
function of the detuning for the case a = .3A and SI = .2. 
The dotted lines show the uncoupled (Q = 0) eigenval- 
ues. The shifts in the asymptotic values of the energy 
levels from the bare trap spectrum {y n — 2n) are due to 
the presence of s-wave collisions. The asymptotic state 
at w —.7 is the bound state of the 'open' channel, which 
is an eigenstate of the trap plus pseudo-potential system. 

In the strong coupling regime, defined as SI 2 ^> 1, 
the width of the resonance is much larger than the trap 
level-spacing, hence many levels can be resonant simul- 
taneously. Thus the low-lying levels all lie very close 
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FIG. 3: Effective scattering length, a e ff, (solid line) and 
molecular fraction, |c n | 2 , (dashed line) as the detuning, 6, is 
swept across resonance. Figures 3a and 3b show the cases 
n = 1 and n = 2, respectively, for the case SI — .2 and a — .3. 
Figures 3c and 3d show n = and n = 1 for SI = 1 and a = 0, 
while figures 3e and 3f show n = and n = 1 for the case 
f2 = 10 and a = 0. The vertical dotted lines mark the location 
of the resonance, while the horizontal dotted lines correspond 
to the analytical result for \c„\ given by Eq. J32V 

to their on-resonance values of v n = 2n — 1. This is 
illustrated in Figure |U which shows the eigenvalue spec- 
trum as a function of detuning for the case il = 10 and 
a = 0. In this regime Eq. I122II is a good estimate for 
the molecular fraction, showing that the molecular am- 
plitude decreases dramatically with increasing coupling 
strength. To understand this effect, we simply make 
the reasonable assumption that in the strong-coupling 
limit all quasi-resonant levels are mixed with equal am- 



plitudes. For il 2 ^> 1, the number of near resonant levels 
is Ni eve i s i=s £1 2 . If we equate the probability for any 
given bare-state to the total probability divided by the 
approximate number of levels we arrive at |c| 2 « l/^ 2 , 
which agrees well with Eq. I|22l) . 

In Figure 3 we plot a e f / and | c„ | 2 versus detuning for 
several cases of interest. In Figs 3a and 3b we show the 
weak-coupling case fl = .2 and a = .3 for levels n = 1 
and n = 2 respectively. The n = 1 case shows a sweep 
(right to left) from the lowest 'unbound' state into the 
bound state in the 'open' channel. The n — 2 case shows 
a transfer from one 'unbound' state to another. As the 
level is swept through resonance we see a broad feature 
in the molecular fraction |c„| 2 , whose maximum value 
is slightly larger than the on-resonance value l|22|l and 
occurs to the right of the resonance. Figures 3c and 3d 
show the intermediate case ft = 1 and a = for levels 
n = and n = 1. We see in the n — 1 case that the 
molecular fraction is significantly reduced compared to 
the weak-coupling regime. Lastly, in Figures 3e and 3f we 
see the strong-coupling case Q — 10 and a = 0, for levels 
n = and n = 1. We see that in the strong coupling 
regime, the scattering length can be tuned from — oo to 
+oo, with a negligible bare- molecular component. 

In conclusion, we see that the effects of pair- 
correlations play a major role in atom-molecule coupling, 
resulting in the appearance of a 1/r singularity in the 
relative wavefunction together with a corresponding de- 
crease in the bare-molecule population. This suggests 
that for molecule formation it is best to have a weak cou- 
pling, while for manipulation of atomic interactions, e.g. 
for BCS pairing of fermions 0, 0], a strong coupling 
will remove the corresponding bare-molecule population. 
In FR the free-space coupling strength is predetermined 
by atomic properties, hence can only be increased by 
decreasing the trap size. In PA, however, the coupling 
strength is readily increased by increasing the laser inten- 
sity. This suggests that laser-induced photoassociation 
may have a significant advantage over Feshbach Reso- 
nance for tuning atom-atom interactions. 
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